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An efficient scheme to compute the geometric entanglement per lattice site for quantum many-body systems 
on a periodic finite-size chain is proposed in the context of a tensor network algorithm based on the matrix 
product state representations. It is systematically tested for three prototypical critical quantum spin chains, 
which belong to the same Ising universality class. The simulation results lend strong support to the previous 
claim [Q.-Q. Shi, R. Orus, J. O. Fja;restad, and H.-Q. Zhou, New J. Phys 12, 025008 (2010); J.-M. Stephan, 
G. Misguich, and F. Alet, Phys. Rev. B 82, 180406R (2010)] that the leading finite-size correction to the 
geometric entanglement per lattice site is universal, with its remarkable connection to the celebrated Affleck- 
Ludwig boundary entropy corresponding to a conformally invariant boundary condition. 
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Introduction. In the last decade, significant progress has 
been made in the investigation of quantum phase transitions 
(QPTs) from a novel perspective-quantum entanglement 1 1]. 
The main idea is to quantify entanglement present in a ground- 
state wave function for a quantum many-body lattice sys- 
tem in terms of a variety of bipartite entanglement measures, 
with some anomaly at a critical point. A remarkable result 
is achieved for the von Neumann entropy that quantifies the 
bipartite entanglement when a finite-size spin chain is parti- 
tioned into two disjoint parts: it scales logarithmically with 
the subsystem's size, with a prefactor proportional to the cen- 
tral charge, a fundamental quantity in conformal field theory, 
as long as the system is at criticality iQ-Iltl- 

Contrary to bipartite entanglement measures, no consensus 
has been achieved regarding multipartite entanglement mea- 
sures, although much effort has been devoted to characterize 
quantum criticality from some multipartite entanglement mea- 
sures. One promising candidate among them is the geometric 
entanglement (GE) [6, 7]. As a holistic measure of the mul- 
tipartite entanglement, the GE quantifies the multipartite en- 
tanglement present in a quantum state wave function. For a 
quantum many-body lattice system, the GE per lattice site is 
shown to be an alternative way to detect critical points E - IJ ] . 

In addition, an intriguing connection between the GE per 
site and the Affleck-Ludwig g factor [12] is established for a 
finite-size spin chain with the periodic boundary conditions at 
criticality. More precisely, the GE per site, e^, scales as 
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with being the lattice size. It was conjectured that the co- 
efficient b in the sub leading term is universal OlSll . This is 
confirmed later in Ref. il4|] for both the transverse Ising and 
XXZ chains at criticality, by relating the coefficient b to the 
Affleck-Ludwig boundary entropy s: 



(2) 



in the sense that the Affleck-Ludwig boundary entropy s ap- 
pears in the GE per lattice site for a periodic chain. However, 
it still remains unclear as to whether or not such a connection 
between the GE per site and the Affleck-Ludwig g factor is 
universally valid, and whether or not there is any criterion to 
judge which g factor, corresponding to a conformally invari- 
ant boundary condition, is chosen for a specific model. There- 
fore, it is desirable to investigate the models belonging to the 
same universality class, to see if the same g factor appears, 
thus providing further insights into the connection between 
the subleading term coefficient b in the GE per lattice site and 
the Affleck-Ludwig g factor. 

One of the main obstacles to address these issues is due 
to the fact that the computation of the GE per lattice site is 
a formidable task for a quantum many-body lattice system, 
because it involves the optimization over all the possible sep- 
arable states, with possibly some constraints arising from the 
translational invariance. However, recent progress in the con- 
text of the tensor network algorithms for quantum many-body 



lattice systems with the periodic boundary conditions Ill5l - ll9ll 



offers us an efficient method to systematically evaluate the GE 
per lattice site for quantum many-body lattice systems. In this 
paper, we first describe how to evaluate the GE per lattice site 
for a finite-size quantum spin chain from a tensor network 
algorithm based on the matrix product state representations. 
Then we present our simulation results for three prototypical 
spin chains belonging to the Ising universality class. 

The geometric entanglement per lattice site. We first in- 
troduce the maximum fidelity A^ax between a quantum pure 
state and all the possible separable and normalized states 
10) of the parties: 



Amax = max K0|iA>l- 

\4>) 
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Here, the boundary entropy s is defined through the Affleck- 
Ludwig g factor: s = In g iS-Sl. This result is surprising. 



The larger A^ax is, the closer a quantum state wave function is 
to a separable state. Therefore, a holistic measure of the multi- 
partite entanglement present in a quantum state wave function 
is defined as 



(4) 
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Since the contribution to Eiij/) from each party is additive, 
£■(1^) scales linearly with A^, for a multipartite system con- 
sisting of parties. Therefore, it is convenient to define the 
GE per party as 
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For a quantum many -body lattice spin system, each lattice site 
constitutes a party. Thus, is the GE per lattice site. Note 
that it is well defined even in the thermodynamic limit. 

A matrix product state algorithm. Let us recall the key steps 
to produce a ground-state wave function from an efficient vari- 
ational algorithm iflTll for a translational-invariant finite-size 
periodic lattice system. First, choose a random state as an ini- 
tial state li/^o). Second, perform the imaginary time evolution 
for If/ro). Thus, we get the evolved state at imaginary time 

T, 



exp(-j/T)|tAo) 
|exp(-i/T)|iAo)l 
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If T —» oo, the ground-state wave function is projected out, as 
long as the initial state is not orthogonal to the real ground 
state. Third, exploit the Trotter- Suzuki decomposition and 
turn the imaginary time evolution operator into a series of two- 
site gates over a time slice 6t, with r = M6t. Therefore, the 
problem becomes how to absorb a two-site gate acted on a ma- 
trix product state. For a periodic spin chain with sites, we 
assume the matrix product state representation is translation- 
invariant under two-site shifts. Thus we only need two three- 
index tensors, Ao and Ae (cf. Fig.[T](i)), to represent a tensor 
network for the ground-state wave function and two one-index 
tensors. Bo and (cf. Fig. [1] (i)), to represent a separable 
state, as shown in Fig. [T](ii) and (iii), respectively. Here, the 
subscripts o and e represent odd and even sites, respectively. 
In this setting, a two-site gate is absorbed by performing a sin- 
gular value decomposition for a matrix contracted from a few 
tensors involving the two-site gate. 

A peculiar feature of the algorithm for a finite-size peri- 
odic system is to take account of the p largest eigenvalues and 
the corresponding eigenvectors of the transfer matrix E(^\^), 
when one evaluates the ground-state energy per site. That is, 
E'!^^^) = Tjf^i Uis'"vi - lif^i M,i"'v,-, where is the i-th largest 
eigenvalue, and v,- indicate the left and right eigenvectors, 
respectively, with D denoting the bond dimension. When the 
energy per site converges, the ground-state wave function is 
generated. Note that the computational cost of the algorithm 
is pD^. 

The geometric entanglement per lattice site from matrix 
product state representations. Now we are ready to compute 
the GE per lattice site. The crucial step is how to maximize 
the fidelity K0|(A)| over all the possible separable states \(p). In 
Fig- [D (iv), we introduce the transfer matrix E(^\^) to repre- 
sent the fidelity between \^) and with -^(^i^f,) constructed 
from two three-index tensors, Ao and A^, and two one-index 
tensors, B* and BJ. Mathematically, the fidelity between a 
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FIG. 1: (color online) (i) Two three-index tensors Aj,", and A'j" are 
attached to odd and even sites respectively, to represent a ground- 
state wave function. Two one-index tensors and 6^, attached to 
odd and even sites respectively, represent a separable state. Here, s 
is the physical index, Z, and r are the inner bond indices, (ii) The 
pictorial representation for a ground-state wave function (iii) 
The tensor network representation for a separable state |(^). (iv) The 
fidelity between the ground-state wave function \tf/) and a separable 
state 10). The transfer matrix E^^i^^ for the fidelity is constructed from 
the tensors A^, A^, B* and S*. (v) The norm for a separable state \(p), 
where the transfer matrix £(010) is constructed from B^, and their 
conjugates. 
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FIG. 2: (color online) (i) The transfer matrix E(^\^) may be approxi- 
mated by the p largest eigenvalues A and the corresponding left and 
right eigenvectors, U and V, of the transfer matrix E^^^). (ii) The 
pictorial representation for the derivative of with respect to S*. 
(iii) The pictorial representation for the derivative of ((p\(p) with re- 
spect to S*. 



ground-state wave function {if/) and a separable state \(p) is ex- 
pressed as: 



f- 



\{m\ 
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For our purpose, we consider F = instead of / itself. To 
maximize F, we compute its logarithmic gradient with respect 
to a one-index tensor B*, 



dlnF 

dB* 



1 dim 1 dim 
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Here, B* is either B* or B*. We take B* as an example to ex- 
plain the updating procedure. The pictorial representation for 
the derivative of In F with respect to B* is shown in Fig. |2] In 
Fig.|2|i), E(^\^) is approximated by the first p largest eigenval- 
ues and the corresponding eigenvectors of E(^\^). That implies 



that E' 



with p specifically depending on 
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m: the larger m is, the less p is. In Fig. |2ii), the derivative 
of with respect to B* is shown. In Fig. |2iii), we rep- 
resent the derivative of {<p\(l>) with respect to B*. Once the 
fidelity gradient is determined, the real and imaginary parts of 
Bo may be updated as follows: 



and 



oBo 



dF 

3Bo = 5Bo + S3(—), 
oBo 



(9a) 



(9b) 



where 6 is the step size in the parameter space, which is tuned 
to be decreasing during the updating procedure. Here, we 
have normalized the tensors B,,, %{dF /dBg) and 3(dF/dBo) 
by setting their respective largest entries to be unity. In exactly 
the same way, we may update the other tensor B^. Repeating 
the procedure until the fidelity per lattice site converges, we 
achieve the closest separable state |0) that ensures its maxi- 
mum fidelity with the ground-state wave function Thus, 
the GE per lattice site follows, as desired. 

The models. We test our scheme by considering three proto- 
typical critical quantum spin chains with the periodic bound- 
ary conditions. Note that the models belong to the same Ising 
universality class. 

The first model is quantum Ising model in a transverse mag- 
netic field on a finite-size ring. The Hamiltonian takes the 
form: 



H 



J'+i] 



(10) 



where cr[^' (a - x, z) are the Pauli spin operators at site /, and A 
is a transverse magnetic field. The model is critical at - I. 

The second model is a spin- 1/2 Ising chain with antisym- 
metric anisotropic and alternative bond interactions. It is de- 
scribed by the Hamiltonian: 

// = 1^ ((1 - (-l)V)^rMcrf " + (afa^^'^ - afaf'^)) , 
(=1 

(11) 

where cr[j^ (a - x,y,z) are the Pauli spin- 1/2 operators at 
site i. The alternative bond interaction is characterized by the 
relative strength r of the exchange coupling, and D, is the z 
direction component of the Dzyaloshinskii-Moriya interaction 
arising from the spin-orbit coupling. The model is critical at 
D- ~ 0.263, if r is chosen to be 0.5 [20]. 

The third model is quantum XYX model in an external 
magnetic field described by the Hamiltonian: 

H^Yj HVI?"" + Aycjfa^^'^ + crl'VI'^" + ha^'^) , (12) 
1=1 

where cr\'} (a - x, y, z) are the Pauli operators at site A,, is 
a parameter describing the rotational anisotropy, and h is an 



external magnetic field. For A,, - 0.25, the critical magnetic 
field is hr ^ 3.206 1I21I42I. 

Simulation results. For the transverse Ising model, we have 
evaluated the GE per lattice site, as shown in Fig. [3a), with 
b - 1.016076, consistent with our previous result in Ref. 



OU. This yields the Affleck-Ludwig g factor g = 0.7032. 
Compared to the exact value ^fixed - 1/ V2 112], the rela- 
tive error is \ig - gfixed)l/§fixed = 1-3 x 10""^. We emphasize 
that, for this model, we have chosen the separable states to 
be translation-invariant under one-site shifts, although there 
is no a priori reason to argue that the GE per lattice site 
does not depend on the unit cell of the separable states un- 
der translational shifts. However, the same g factor is yielded, 
even if the separable states are translation-invariant under two- 
site shifts. In fact, we have evaluated the GE per site with- 
out any translation-invariant assumption, the same g factor is 
achieved. Thus, we conclude that only the smaller g factor, 
corresponding to the stable conformally invariant boundary 
condition |24], i.e., fixed boundary condition, is involved. 

In Fig. Ob), we plot the GE per lattice site, &n, as a 
function of the number of the lattice sites, A^, for the spin- 
1/2 Ising chain with the coupling parameters r - 0.5 and 
- 0.263. The finite-size scaling for the GE per lattice 
site follows Sn - a + b/N + f/N^, where the coefficients are 
a = 0.083320, b = 0.966353 and / = -1.398155. It yields 
the Affleck-Ludwig g factor g - 0.715330. Compared with 
the exact value ggxed - 1 / V2 lfl2ll for the conformally in- 
variant fixed boundary condition, the relative fitting error is 
l(^ - gfixed)l/gfixed = 1.1 X 10"^. For this model, we have cho- 
sen the separable states to be translation-invariant under two- 
site shifts, given the alternative bond interaction. However, 
we have evaluated the GE per site without any translation- 
invariant assumption, the same g factor is achieved. Thus, 
we conclude that only the smaller g factor, corresponding to 
the stable conformally invariant boundary condition |24J, i.e., 
fixed boundary condition, is involved. 

In Fig. [3|c), we plot the GE per lattice site, Sn, as a func- 
tion of the number of the lattice sites, A^, for quantum XYX 
model in an external field, with A,, = 0.25, where the size 
A^ is chosen from 10 to 100. Remarkably, it splits into two 
branches for a large enough chain size: one is for the sepa- 
rable states that are translation-invariant under one-site shifts, 
and the other is for the separable states that are translation- 
invariant under two-site shifts. The exact diagonalization for 
small-size chains up to 24 is also performed to evaluate the GE 
per site without any translation-invariant assumption. How- 
ever, it only reproduces the second branch, as seen in Fig.[3jc). 
A proper explanation for this unexpected result is that, for this 
model, the translational invariance under one-site shifts con- 
stitutes a true constraint to separable states. That is, the clos- 
est separable state with the one-site unit cell is different from 
the closest separable state with the two-site unit cell. Thus 
the maximum fidelity Amax is larger for the latter, implying 
that the GE per lattice site is smaller, as long as the chain 
size is large enough. This explains the presence of a thresh- 
old in the chain size: the GE per site is the same for both 
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FIG. 3: (color online) Main: A scaling relation between the GE per 
lattice site Gjy and the chain size N, for quantum Ising model (a), 
with the size A' ranging from 8 to 120, for quantum Ising model with 
antisymmetric anisotropic and alternative bond interactions (b), with 
the size N ranging from 8 to 150, and for quantum XYX model (c), 
with the size N ranging from 10 to 100. The data are fitted into fi^. = 
a+b/N+f/N-, witha = 0.030960, b = 1.016076 and / = -1.712521 
in (a) and a = 0.083320, b = 0.966353 and / = -1.398155 in (b). 
In (c), we have a = 0.056892, i = 0.001581, and / = -0.879631, 
and a = 0.010193, i = 1.009600 and / = -5.005750 , respectively, 
for the upper and lower branches. Here, the upper (lower) branch is 
for the separable states that are translation-invariant under one-site 
(two-site) shifts, if the chain size is lager than a threshold. Indeed, 
the GE per site is the same for two cases, if the chain size is less 
than the threshold. In addition, the exact diagonalization results for 
the GE per lattice site are displayed in five-pointed stars, up to the 
chain size 24. The data match very well. Inset: The relative fitting 

error E^r = (e^"'" - ^n)/^n ^''^ ^^^^ ^-^ ^ ^^^^ ^^^^ 

1 .8 X 10"^ in (b), and less than 1 .5 x 10"'' in (c), where is the value 
extracted from the fit and fij"" is our simulation value for each A'. 



the one-site and two-site unit cells, if the size is less than the 
threshold. Two sets of the data are separately fitted into the 
scaling function &ff = a + b/N + f/N^, with the coefficients 
a = 0.056892,/? = 0.001581, and / = -0.879631, and a = 
0.010193,^7 = 1.009600 and / = -5.005750, respectively, 
yielding gxYX = 0.9994 and gxYX = 0.7048. They are con- 
sistent with the exact values gfree — 1 and ggxed = 1 / V2 
for conformally invariant free and fixed boundary conditions 
in the Ising universality class. 

Summary. In this paper, we have developed a scheme to 
efficiently compute the GE per lattice site for quantum many- 
body spin systems on a periodic finite-size chain in the con- 
text of a tensor network algorithm based on the matrix product 
state representations. The computational cost does not depend 
on the chain size. A systematic test is performed for three pro- 
totypical critical quantum spin chains belonging to the same 
Ising universality class. The simulation results lend strong 
support to the previous claim that the leading finite-size cor- 
rection to the GE per lattice site is universal 1 13], with its re- 
markable connection to the celebrated Affieck-Ludwig bound- 
ary entropy corresponding to a conformally invariant bound- 
ary condition 1 14]. For all the models tested, the simulated g 
is compared to the exact g factor from conformal field theory, 
with the relative error less than 1.1 x 10"^. It appears that the 
boundary entropy corresponding to the smallest g factor is al- 
ways involved. This is somewhat expected, since this g factor 
characterizes a stable fixed point along a boundary renormal- 
ization group flow, according to the Affleck-Ludwig g theo- 
rem 1 12]. Remarkably, for quantum XYX model in an exter- 
nal field, either conformally invariant free or fixed boundary 
condition appears, depending on the one- or two-site unit cell 
of the translation-invariant separable states. 
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